Abstract. We discuss, both analytically and numerically, 
the behavior of magnetic pressure and density fluctuations 
in strongly turbulent isothermal magnetohydrodynamic 
(MHD) flows in "1+2/3" dimensions, or "slab" geome- 
try. We first consider "simple" MHD waves, which are the 
nonlinear analogue of regular MHD waves, and have the 
same three modes, slow, fast and Alfven. These allow us 
to write equations for the magnetic field strength i? as a 
function of density p for the slow and fast modes, show- 
- - - ing that the two have different asymptotic dependences 
' of the magnetic pressure vs. p. For the slow mode, 
i?^ ~ ci — C2P, while for the fast mode, ~ p^ . We also 
\ perform a perturbative analysis to investigate Alfven wave 
^jpj' pressure, recovering the results of McKee and Zweibel that 
^ B^ ~ p2, with 7e ~ 2 at large Ma, 7e — 3/2 at moder- 
■ ate Ma and long wavelengths, and 7e — 1/2 at low Ma- 
This wide variety of behaviors implies that a single poly- 
tropic description of magnetic pressure is not possible in 
I I general, since the relation between magnetic pressure B^ 
^ ■ and density is not unique, but instead depends on which 
C<1 I mode dominates the density fluctuation production. This 
. in turn depends on the angle 9 between the magnetic field 
"J I and the direction of wave propagation and on the Alfvenic 
. Mach number Ma- Typically, at small Ma, the slow mode 
' dominates, and B is anticorrelated with p. At large Ma, 
both modes contribute to density fluctuation production, 
and the magnetic pressure decorrelates from density, ex- 
O ^ hibiting a large scatter, which however decreases towards 
I higher densities. In this case, the magnetic "pressure" does 
^ not act as a restoring force, but rather as a random forc- 
' ing. These results have implications on the probability 
Ci • density function (PDF) of mass density. The unsystem- 
' atic behavior of the magnetic pressure causes the PDF 
to maintain the lognormal shape corresponding to non- 
magnetic isothermal turbulence, except in cases when the 
^ slow mode dominates, in which the PDF develops an ex- 
" " ' cess at low densities because the magnetic "random forc- 
ing" becomes density dependent. Our results are consis- 
tent with the low values and apparent lack of correlation 
between the magnetic field strength and density in sur- 
veys of the lower-density molecular gas, and also with the 
recorrelation apparently seen at higher densities, if the 
Alfvenic Mach number is relatively large there. 
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1. Introduction 

The cold molecular gas in our Galaxy is generally believed 
to be magnetized and turbulent (see, e.g., the reviews by 
Dickman 1985; Scalo 1987; Heiles et al. 1993; Vazquez- 
Semadeni et al. 2000). However, at present, the actual 
strength of the magnetic field in molecular clouds and, as 
a consequence, its dynamical importance, are still a mat- 
ter of debate, with opinions ranging from considering it 
crucial (e.g., Crutcher 1999) to moderate (e.g., Padoan & 
Nordlund 1999), although in general, the observational ev- 
idence so far appears inconclusive (e.g., Bourke et al. 2001; 
Crutcher, Heiles & Troland 2002). Therefore, it is impor- 
tant to continue gathering information both observation- 
ally and theoretically on the distribution and strength of 
the magnetic field in the ISM and, in particular, molecular 
clouds. 

In general, magnetic fields and turbulence are believed 
to provide support against the self-gravity of the clouds, 
which typically have masses much larger than their ther- 
mal Jeans mass. In this context, the pressure provided by 
magneto- hydrodynamic (MHD) motions in the clouds has 
been investigated in a variety of scenarios. Since obser- 
vations often suggest that there is near equipartition be- 
tween turbulent and magnetic energies in molecular clouds 
(e.g., Myers & Goodman 1988; Crutcher 1999), the mo- 
tions are often considered to be MHD waves (Arons & Max 
1975). Until recently, Alfven waves were favored, as these 
were expected to be less dissipative than their slow and 
fast counterparts, but recent numerical simulations (Mac 
Low et al. 1998; Stone, Ostriker & Gammie 1998; Padoan 
& Nordlund 1999) have shown that supersonic MHD tur- 
bulence decays almost as rapidly as purely hydrodynamic 
turbulence because of strong coupling between the MHD 
modes, and so the rationale for the preponderance of the 
Alfven mode appears less clear now (but see Cho, Lazar- 
ian & Vishniac 2002 for an argument against fast decay). 
It is thus necessary to investigate all MHD modes as a 
source of pressure in a general context in fully turbulent 
regimes. 



The functional form that the magnetic pressure may 
take as a function of density is also important for the 

statistics of the density fluctuations, as it was shown by 
Passot & Vazquez-Semadeni (1998; see also Nordlund & 
Padoan 1999) that for flows in which the pressure P ex- 
hibits an effective polytropic behavior, such that it scales 
with density p as P ^ p^", with 7e the effective poly- 
tropic exponent, the shape of the density probability den- 
sity function (PDF) depends on the specific value of 7e. 
The PDF takes a lognormal form for 7e = 1, but it devel- 
ops an asymptotic power-law tail at low (high) densities 
for 7e > 1 (7e < !)• Thus, the value of 7e produced by the 
various sources of pressure is important for the density 
statistics and the overall dynamics of the flow. 

In an analysis of the pressure produec;(l by Alfven 
waves, McKee & Zweibel (1995, hereafter MZ95) con- 
cluded, based on an analysis by Dewar (1970), that it is 
isotropic, and that, for flows undergoing slow compres- 
sions, 7o = 3/2, while for strong shocks, 7c = 2. Vazquez- 
Semadeni, Canto & Lizano (1998) studied the evolution 
of the velocity dispersion a in gravitationally collapsing 
flows, as a means of determining the density dependence of 
the "turbulent pressure" Pturb, defined in such a way that 

= dPtuih/dp. They found that the for slowly collaps- 
ing magnetized cases, 7e ~ 3/2, while for rapid collapse, 

7e ~ 2. 

More recent numerical work has directly plotted the 

magnetic strength B or the magnetic pressure in 
three-dimensional numerical simulations of both isother- 
mal (Gammie & Ostriker 1996; Padoan & Nordlund 1999; 
Ostriker, Stone & Gammie 2001) and non-isothermal flows 
(Kim, Balsara & Mac Low 2001; MacLow et al. 2001). 
The isothermal simulations should be reasonable models 
of molecular clouds, while the non- isothermal ones are cast 
as models of the ISM at large. Gammie & Ostriker (1996) 
report an almost constant (on average) magnetic field 
strength as a function of density in non-self-gravitating 
runs. Padoan & Nordlund's (1999) plot of i? vs. p shows a 
scatter of over two dex in the low density range for large 
values of the Alfvenic Mach number, with the scatter de- 
creasing towards larger densities. On the other hand, there 
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is less scatter at low Alfvenic Mach numbers (Ma), but it 
is more uniform throughout the density range. Those au- 
thors note that the upper envelope of their B distributions 
closely matches the observational scaling B ~ p^/^ found 
in cases when the magnetic field is detected (e.g., Crutcher 
1999; Crutcher et al. 2002), although they remark that in 
a large number of cases only upper limits or complete non- 
detections are obtained. Such scaling also corresponds to 
the high-S tail of the magnetic strength distribution in the 
high-A/a simulation by Padoan & Nordlund. It is worth- 
while to additionally note that the bulk of the points in 
their B distribution does not follow a simple power-law 
scaling with p, but instead appears to curve up, being 
first nearly constant with p and then starting to rise. Os- 
triker et al. (2001) find similar results, except that these 
authors only attempt to fit the high-density part of the 
B distribution to observations, finding slopes between 0.3 
and 0.5, which again are not too different from the ob- 
served scaling. 

In non-isothermal simulations, Mac Low et al. (2001) 
have plotted the magnetic pressure vs. the thermal pres- 
sure, finding again a large scatter of up to 6 orders of mag- 
nitude in their supernova-driven simulations, while Kim et 
al. (2001) show plots of B vs. p, finding slopes of ^ 0.4, al- 
though with scatter of one order of magnitude at high den- 
sity, and two orders of magnitude at lower densities, the 
plot actually appearing more like two segments with very 
different slopes. Finally, in simulations of thermal conden- 
sation triggered by strong compressions in the presence of 
a uniform magnetic field, which are in several aspects sim- 
ilar to gravitational collapse, Hennebelle & Perault (2000) 
found that the magnetic field does not necessarily increase 
together with the density. 

Recent theoretical work (e.g., Lithwick & Goldreich 
2001; Cho et al. 2002) has not specifically addressed the 
issue of the magnetic pressure-density correlation, as it 
has focused mainly on the spectral properties of moder- 
ately compressible MHD turbulence, as a consequence of 
mode coupling. Nevertheless, Cho et al. (2002) mention in 
passing that their results lead them to not expect a signif- 
icant correlation between the magnetic pressure and the 
density. 

In an attempt to understand the physics underlying 
the above experimental and observational results, in this 
paper we study the dependence of magnetic pressure with 
density as a consequence of nonlinear MHD wave propa- 
gation in isothermal MHD turbulence in "1-1-2/3" dimen- 
sions, also referred to as a "slab" geometry. This choice 
allows us to isolate the relative importance and effects 
of the relative orientation of the magnetic field and wave 
propagation directions. We first consider the problem ana- 
lytically, finding asymptotic relations between the density, 
magnetic field and velocity field fluctuations, which show 
the effectiveness of the slow and fast modes as sources of 
density fluctuations under a variety of conditions (S^). We 
then discuss Alfven wave pressure using a perturbation 



analysis, comparing with the results of MZ95 (§H). These 
results are then tested by means of numerical simulations 
in the same geometry (§^), which support the scenario 
derived from the analysis of the equations. Finally, in ^ 
we present a summary and some discussion of our results 



(§5.1), including implications for the interpretation of ob- 
servational results (§^.2[). 



2. Properties of simple waves 

The equations governing, in the magnetohydrodynamic 
(MHD) limit, the onc-dimcnsional motions of a plasma 
permeated by a uniform magnetic field Bq in a slab ge- 
ometry read 



dp _^ djpu) 
dt dx 



du 
'dt 



du 





1 d 



— +u— = —— \-^ + ^]+f, 



dx 



\b\' 



dv dv 

L u = 

dt dx M^p dx 
db 9 , ^ dv 



pdx 2M2 
bx db 



f 



(1) 
(2) 
(3) 

(4) 

where x is the direction of propagation, nondimensional- 
ized by a typical length L. The velocity components u 
(along the x axis) and v = Vy + iv^ (along the y and z 
axes) are normalized by a velocity unit Uq, the mass den- 
sity p by a reference density pQ, and the magnetic field 
components — cos 9 and b = by + ib^ by Bq . Note that 
bx is constant and that bz{t = 0) = sin 6*, where 9 is the 
angle between the direction of propagation and that of 
the initially unperturbed ambient magnetic field. Time t 
is measured in units of L/Uq. An isothermal equation of 
state is assumed and we denote by Cg the constant sound 
speed. Two non-dimensional numbers can be defined, the 
sonic Mach number Mg = Uq/cs and the Alfvenic Mach 
number Ma — Uo/va, where Va = i3o/(47r/9o)^^^ is the 
Alfven speed of the unperturbed system. The plasma beta 

is here defined by /3 = . The above system of equations 
also contains a driving in the form of random accelera- 
tions fx (acting on the component u of the velocity) and 
f = fy + ifz, (acting on the components v). 

In this section we shall derive some properties of the 
so-called simple MHD waves (see, e.g.. Landau & Lifshitz 
1987, §101) and thus we assume fx = f = 0. In the case 
where the basic state is perturbed by infinitesimal dis- 
turbances, the solutions are superpositions of linear plane 
travelling waves. These plane waves are monochromatic 
and their profile does not change in time, with all quanti- 
ties only depending on the combination x±Ct. The prop- 
agation velocity C is a constant that identifies with one of 
the three possible roots of the dispersion relation, namely 
the slow, Alfven or fast velocity. For a particular plane 
wave solution each perturbed quantity can be expressed 
as a function of a chosen one, for example the density p. 



4 



Passot & Vazquez-Semadeni: Density-magnetic pressure correlation 



In the case of finite-size perturbations these relations 
do not hold, but it is nevertheless possible to look for par- 
ticular solutions that have the property that all quantities 
are only functions of any single one of them, as in the 
linear theory. These particular solutions, which general- 
ize the linear plane wave solutions to the case of finite 
disturbances, are called simple waves. 

Following Landau and Lifshitz (1987, §105; see also 
Mann 1995), we recall here how to derive the relevant 
equations that determine simple wave profiles. The prin- 
ciple is illustrated on the equation for mass conservation, 

that can be rewritten in the form 
d{p,x) d{pu,t) 



0, 



(5) 



d{t,x) d{t,x) 
after transforming the partial derivatives into Jacobians. If 
one now chooses (t, p) as new independent variables, one 
has to multiply the above equation by . Assuming 
that all dependent variables only depend on p, and writing 
U = ^ one gets, after expanding the Jacobians, 

-U + p^+u^O, (6) 
dp 

or, more generally, 

- Vdp + pdu = 0, (7) 
where V = U—u denotes the wave speed. Each field is thus 
a function F{x — {u + V)t) where u and V are functions 
of e.g. p. Each point of the wave is traveling with its own 
velocity, leaving the possibility for wave steepening and 
the subsequent formation of a discontinuity. 

Using the same procedure, equations read 
1 1 \b? 



-Vpdv 



rdb = 



(9) 



M? 

h 

-Vdb + bdu - b^dv = 0. (10) 
The system (0)-(|l^) has non-trivial solutions if the 
wave speed V is given by 

1 



(B^+Pp) 1± 



if3blp 



(11) 



^/3p)V 

''"^ Recall that 



2M2p ^ \^ V (52 

or if it equals the Alfven speed Va — ijj—^^ 
we denote by ^ b^ + |6p the total magnetic intensity. 

The latter root is associated with the circularly po- 
larized Alfven simple wave, which is non-compressive and 
has |6p = est. The solutions ( pi] ) are associated with fast 
and slow simple waves. The speeds of the linear fast and 
slow magnetosonic waves are recovered by taking p—1 and 
b = i sinO. 

We are thus led to the following system of equations 
for the wave profiles 

dp p 

d d B^ 



dp 2 dp 2 ^ ° 



P). (13) 

Equation (|13| ) can also be rewritten as dP/ dp = V^, where 
P denotes the total pressure -I- j^. Equations ( pi] ) 



and (|12[)-(p3[) can be solved numerically but it is advan- 
tageous to search for asymptotic solutions in some limits. 
First of all, when f3 ~ 0, i.e. in the absence of thermal 
pressure, it is found that V- = (the slow wave does 
not propagate) and = -^r^, which is the Alfven speed 
based on the total magnetic field intensity. When bx — 0, 
i.e. for a propagation perpendicular to the ambient mag- 
netic field, we again have V- = and now = -^jgs^- 
In both limits, for the slow waves P = const, while for 
the fast waves B^ (x p^. These relations are in fact more 
general as will be seen below. 

In the case where 4/36^p <C {B^ + (ipY one obtains 
bl 



M^{B^ 
B^ 

MIp + M| 



/3p) 



(14) 



(15) 



The above assumption only fails when ^ = ^ is of or- 

der unity together with rj ~ -^4- small, i.e., when bx is 
not too small, for j3p of order unity and small field dis- 
tortions. As an example, the above approximations fail 
for sound waves propagating along the magnetic field, for 
which arbitrary density pertubations can develop on an 
unperturbed magnetic field. 

Equation (|l3|) together with the approximation (15) 
leads to B^ oc p^ for the fast wave. Using eq. ( p^ ) for the 
slow speed, one is led to the following ordinary differential 
equation 



whose implicit solution reads 

e + | + ln|l-e-r?|-C', (17) 

C being an arbitrary constant. For large enough ^ and 
?/, the solution can be approximated by 77 = C" — 2^. 
In physical terms it can be rewritten as P const. For 
smaller values of ^ and r\ it can be shown that P is an 
increasing function of p, with slope 1/2 when r/(^ = 0) = 1, 
and even larger slope at smaller values of ^ and r\. 

In order to interpret the numerical simulations of Sec- 
tion 4, it is also useful to investigate under which con- 
ditions do the slow or the fast waves dominate the pro- 
duction of density fluctuations. A criterion that can be 
used is the relation between the total velocity vector 
U = (u,Vy,Vz) of the perturbation (in some way related 
to the fluid displacement) and the density perturbation. 
From eqs. ([7])-(pl]), it follows that, for small enough per- 
turbations, the density fluctuation Ap and the velocity 
fluctuation AU± = (Au±,Avy±,Avz±) are related by 
Ap 



Au± 

and 
Av± 



V±- 



bxb 



:{VlMl-(i) 



Ap 



These relations are not valid 
Alfven waves where |6| is constant. Using eq. 



P 

for 



circularly polarized 
~T|), it is 
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easy to verify that AU+ 



i.e. when by = 0, = sin 



±Bo 

II Bo 



AU_ = 0. In the hnear case, 
' and (0=1, one finds that, in 

the hmit (3^0 

(X sin 6x — cos 9z 
AU_ oc cos 6x + sin 9z 
while in the limit (3 ^ oo 
AU+ oc X II k 
AU_ cxz ±k 
where k is the wavevector of the linear wave. 

Denoting X± — V^M^ and using the dispersion rela- 
tion pXl - {B^ + (ip)X± + [ibl = 0, one finds 

-l^[B^X^+(3pX^-2pbl) 

In the limit /3 ^ one has X 
that 



|AC/±| 



(18) 



= -2- and X_ = 



so 



|AU+p = 



Ap 



Mi \ p ) 



while in the limit /3 cxd, 
P_(Ap_^'' 
Ml \ p , 



= /3, X_ 



= and 



lAU+M = 



|AU_M = 



M2|6P ^ p ) 



From these relations one can conclude that at small /3, 
density fluctuations are mostly created by the slow mode 
(since small values of |AU_| are sufficient to obtain ^ ^ 
1). Conversely, the fast mode can more easily generate 
density fluctuations at large [3. For intermediate values of 
/3, conclusions may depend on the magnitude of b^. For ex- 
ample, at parallel propagation and large field distortions, 
the slow mode tends to be the most efficient to produce 
density fluctuations, except at large density. 

3. Alfven wave pressure 

In the previous Section we considered waves propagating 
on a uniform background. In a turbulent medium, the sit- 
uation is usually more complex, as all types of waves are 
mixed. Another case that still remains simple enough to 
be analytically tractable, consists in studying the prop- 
erties of MHD waves propagating on top of a circularly 
polarized, parallel-propagating Alfven wave. These Alfven 
waves are exact solutions of the MHD equations (Ferraro, 
1955) and can be taken of arbitrary amplitude. They read 
6o = exp [—ia{kQX — toot)] with pa — 1, bxo = 1, uo = 
and vo — —bo/Ma- The polarization of the wave is deter- 
mined by the parameter cr, with (t = +1((t = — 1) for a 
right-handed (left-handed) wave. The dispersion relation 
reads loq = ko/Ma- 

Let us now consider perturbations of the form 

^ = -(_ fj^^-iTlik+ko)x-{uj+uio)t] 

-|-6_e'^*'^[^'^~'^°^'^~^"~"°-'*l (19) 



p=l + p'e-»'^('=^-'^*) + . 
u = u' e-"'^^'-'-'^'^ + c.c. 



4_y ^+i(y[(k-ka)x-(Lj-LJo)i\ 

The linearized equations read 
Lop' — ku' = 
I k I k 



(w + UJo)v+ 
(w " UJo)v^ 



^ 2MI 



kobp 
Ma ' 
koba 



{bob*_ 



bib. 



(20) 
(21) 

(22) 

(23) 
(24) 

(25) 

(26) 

(27) 
(28) 

After some algebra, we obtain the following dispersion re- 
lation 



-P = Jp{k + kQ)b, 



Ma 

uJobo 



Ma Ma ' 

{uj + uJo)b+ — (fc + fco)^ow' = — (fc + fco)w+ 
{lo — LUo)b^ — (fc — k{))b{)U * = — (fc — ka)v^ 



Ml 



Ml-^^ 



Ml 



Ma) \ Ma Ml Ml) ^ ' 

Three different limiting cases can be considered and easily 
identified after rewriting the dispersion relation as 

^^^ = ^77 + 
Ml 



ffBl 
Ml 

= 0. 



fc2S2 I ^ 



k 2 



fc2 



M3 



Ml 



(30) 



In the case Ma oo, i.e. in the case of a weak back- 
ground magnetic field, magnetic tension is negligible com- 
pared to field stretching and the dispersion relation ap- 



proximates to UJ 



k^ 
17? 



M2 



If one assumes a polytropic 



dependence of the magnetic pressure on the density in the 
form |6p/2 cx p"^ and linearizes eqs. (|l])-(||), a direct com- 
parison of the resulting dispersion relation with the above 
approximation leads to 7 = 2. This is the case, mentioned 
by MZ95, in which the perturbation is very rapid com- 
pared to the speed of the Alfven wave. This behavior of 
the magnetic pressure is easily analyzed directly from the 
orginal MHD equations. Indeed, in the limit where Ma is 
very large, the term on the right-hand-side of eq. (^) can 
be dropped and this equation leads to 



d\b\' 



d 



d 



dt dx 
Together with eq. 



dx 

one obtains 



(31) 
(32) 



confirming the result 7 = 2 in this case. 

In the case Ma 0, i.e. in the case of a strong back- 
ground magnetic field, magnetic tension is dominant and 
the Alfven wave is very rapid. When « and fc ^ fcp, 
i.e. when the perturbation is very slow, or "quasi-static" , 
the dispersion relation approximates to lu'^ 



k^ 



k'B^ 
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giving 7e = 
proximation. 

Finally, for Ma finite but 
quasi- uniform perturbations, 

4M2 



also derived by MZ95 using a WKB ap- 



proximates to cj « -|W 



fc « and Lo ajQ, i.e. for 
the dispersion relation ap- 
giving 7e = i , recovering 



the prediction of MZ95 for this case 



Whereas the predictions based on the WKB approach 
are recovered in this purely one-dimensional linear anal- 
ysis, it should be noted that they are probably not valid 



in a fully three-dimensional situation. For example, MZ95 



argued that the Alfven wave pressure is isotropic. A linear 
stability analysis analogous to the previous one but per- 
formed in the long-wave limit and for perturbations prop- 
agating exactly perpendicular to the background Alfven 
wave shows that the Alfven wave exerts a pressure differ- 
ent from the ID case, thus invalidating the isotropy result. 
This pressure can even be negative and one recovers the 
polytropic index 7e = 3/2 only in the limit of large am- 
plitude background wave and for /3 = (Passot & Gazol, 
in preparation). 

4. Numerical simulations 

In order to test the range of validity of the analytical 
predictions, we have performed numerical simulations of 
eqs. (0)-(^ in a periodic domain using a pseudo-spectral 
method based on Fourier expansions. In order to numer- 
ically handle the formation of strong shocks, dissipative 
terms of the usual form, namely — i-^ and — -i-r, are 
added to the right-hand-side of eqs. (||) and (^, respec- 
tively, and a magnetic diffusion to the right-hand- 
side of eq. (0). In addition, it was found necessary to 
add a mass diffusion in the form A'-rg;^ to the right- 
hand-side of eq. (^. This term preserves mass conserva- 
tion and allows handling of strong shocks. If kept small 
enough, it does not significantly modify the dynamics and 
in particular it does not alter the statistical conclusions 
we shall present. Except for one set of simulations dis- 
cussed below, the forcing (actually an acceleration) act- 
ing on these equations is applied on Fourier modes 1- 
19, peaked at wavenumber 8 with amplitudes propor- 
tional to fc"* exp(— /c^/32) and Gaussian-distributed ran- 
dom phases. A forcing is also applied on mode in order 
to ensure momentum conservation. The random phases 
are changed every 0.003 time units. A state of constant 
density and zero velocity and magnetic field fiuctuations 
is taken as initial conditions. A resolution of 4096 grid 
points is used and typical values for the diffusion coeffi- 
cients are fJ.i = iJ,t = V = 0.003 and fj,r = 6. x 10~^. In the 
simulations, three main parameters are varied, namely the 
angle 9 and the Mach numbers Mg and Ma, but only 6 and 
Ma are actually important as we only consider the high- 
Mg limit. For each simulation, we compute r.m.s. values 
of the sonic Mach number Mg = Mg {U^} , where the 
brackets (resp. bars) denote time (resp. spatial) averaging. 



Similarly, the r.m.s. Alfvenic Mach number and the effec- 

1/2 

five beta of the flow are defined as Ma = Ma ' 



and (3 — 



Ml 



'W 



Finally, the density fluctuations and 

1/2 

and 



field distortions are defined as ^ = ^(p — 1)" 
^ = ((62 + (6^-sin6')2)2^^^^ respectively. 




Fig. 1. Log- log scatter plot of the magnetic pressure ver- 
sus density for a simulation with ~ 7r/2, = 4.12, 

Ma = 5.16, /3 = 1.52, ^ = 0.623, ^ = 0.618 and 
(7 — 0.054. The line segment has a slope equal to 2. 



The simplest case corresponds to perpendicular prop- 
agation [9 = 7r/2) where only fast waves can propagate. 
In that case, we expect magnetic pressure to be very well 
correlated with density, with a dependance of the form 
5p oc /5^, whatever the value of Ma- This prediction is 
confirmed by the simulations as examplified in Fig. 1, 
which shows the logarithm of the two-dimensional his- 
togram /i [ln(/9), ln(|&p/2M^)] of the magnetic pressure 
versus density for a run with Ma = 5.16, adding points 
from temporal snapshots taken every 0.1 time units from 
t — 100 to t = 10''. More precisely we display in grey 
scale (with black (resp. white) denoting the smallest (resp. 
largest) value) the logarithm of the number of points (in 
the space-time sample) in a given interval of |6|V2M2 and 
/5, plotting only the points where the histogram is larger 
than 90% of the maximum at fixed p. We also calculate 
an average histogram dispersion in the following way. For 
each density p, we define the mean of the magnetic pres- 
sure logarithm as 



Pniip) 



vh [ln(p), x] dx/H 



(33) 
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and the dispersion as 
a{p) = 



X - Vra{p)f h [\n{p),x] dx/H] 



1/2 



(34) 



where H = h{p, x)dx. An average dispersion a is eval- 
uated by taking the mean of &{p) for \n{p) varying by 10% 
about its central value. 





Fig. 2. Log- log scatter plot of the magnetic pressure ver- 
sus density for two simulations with cosO = 0.1. The pa- 
rameters of the left panel, corresponding to a case with 
forcing on all three velocity components are Ms = 6.63, 

Ma = 0.478, /3 = 0.00685, ^ = 2.69, ^ = 0.323 and 
a — 0.18. For the right panel, where the forcing is ap- 
plied on the V velocity components only, Ms ~ 7.27, 



Ma = 0.527, 
a = 0.059. 



0.00728, ^ = 



3.55, 



^ = 0.324 and 



As the angle 9 is decreased from 7r/2, the behavior is 
very sensitive to M^. When Ma <^ 1 the field lines are 
still almost unperturbed but the density fluctuations are 
very different from the case 9 = 71/2 because for small 
(3 the latter are mostly created by slow waves. In that 
case again we expect the magnetic pressure to exhibit lit- 
tle scatter, as the field is only slightly perturbed, and to 
be roughly anti-correlated with density, as the total pres- 
sure remains roughly constant. This is indeed the case, 
as seen in Fig. 2a. where the distribution of |&|^/2M^ 
and p has a small scatter. In addition, a clear anticor- 
relation is observed at high density when no forcing is 
applied on the a;-component of the velocity (Fig. 2b). 
In the case of purely perpendicular propagation the den- 
sity structures are those typically observed in neutral tur- 
bulence with a polytropic index 7 = 2, i.e. rather flat- 
topped structures (plateaux) separated by shocks that 
keep interacting, with a preeminence of low density re- 
gions (Passot & Vazquez-Semadeni 199S). In contrast, for 



nearly perpendicular propagation and Ma <C 1, the den- 
sity structures are composed of large peaks that oscillate 
while avoiding collision (Fig. 3a). The probability density 
function (PDF) of the density shown in Fig. 3b shows 
an extended tail at large density with a significant ex- 
cess at small density. Since the field strength is roughly 
constant and independent of p, the magnetic pressure 
term in eq. (H) acts like an random acceleration whose 



4000 




Fig. 3. Top: Snapshot of logj^g P (solid) and logj^Q |6p (dot- 
ted) for the run of Fig. 2a. Bottom: Density PDF for the 
run of Fig. 2a in log-log coordinates. The dashed line 
shows a fit of the curve with a parabola corresponding, 
in these coordinates, to a lognormal distribution. 



strength increases as p decreases since |6p is roughly con- 
stant at small density (see Fig. 2a). Whereas the density 
field of an isothermal gas stirred by a random accelera- 
tion has a lognormal distribution, when the acceleration 
depends on the density, the PDF ceases to be lognormal 



(Passot & Vazquez-Semadeni 1998). The lift up of the left 



tail seen in Fig. 3b is thus expected in the present case. 

When the angle is still slightly smaller than 7r/2 but 
with Ma » 1, the magnetic field lines undergo large fluc- 
tuations and both slow and fast waves contribute to the 
production of density fluctuations. As seen in Fig. 4, and 
Fig. 5a, magnetic pressure and density are now positively 
correlated (magnetic pressure roughly follows a polytropic 
law with 7 = 2, corresponding to the fast mode), but the 
dispersion is larger than in the case of a small Alfvenic 
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Fig. 4. Log-log scatter plot of the magnetic pressure ver- 
sus density for a simulation with forcing on the three ve- 
locity components, cos^ = 0.1, Mg = 5.73, Ma = 7.29, 

(3 = 1.53, ^ = 1.32, ^ = 2.17 and a = 0.61. The line 
segment has a slope equal to 2. 



Mach number, indicative of the additional contribiition 
of the slow mode. As a result, the stirring due to mag- 
netic pressure is statistically non-correlated with the den- 
sity and the density PDF is expected to be much closer to 
a lognormal as seen in Fig. 5. 

When the angle is intermediate between parallel 
and almost perpendicular propagation, the distinction be- 
tween the small and the large Alfvenic Mach number cases 
is not as clear. As shown in Figs. 6a-b for runs with an 
angle 9 « 7r/4 and for Ma = 0.521 and Ma = 2.48 respec- 
tively, there is a larger scatter of the points for small to 
intermediate values of the density at large Alfvenic Mach 
number. In this case a positive correlation between mag- 
netic pressure and density is noticeable at large density 
(Fig. 6b). 

For parallel propagation with perpendicular forcing, 
the behavior is again strongly dependent on the Alfvenic 
Mach number. When Ma » 1, the magnetic field is 
strongly distorted. The forcing excites slow modes which 
develop into shocks, inside which large density clumps 
form. These clumps, located at local minima of the mag- 
netic field intensity, are separated by regions of large mag- 
netic pressure (Fig. 8a). As a result, they can hardly ap- 
proach each other but rather oscillate about their mean 
position, at a frequency close to that of the fast wave of 
the local total magnetic field. The signature of the slow 
wave dominance is the rather small scatter and the anti- 
correlation between magnetic pressure and density in com- 
pressed regions, as seen in Fig. 7a corresponding to a case 



2 
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Fig. 5. Top: Snapshot of log;^Q p (solid) and log^g |6p (dot- 
ted) for the run of Fig. 4. Bottom: Density PDF for the 
run of Fig. 4 in log-log coordinates. The dashed shows a 
lognormal fit. 



with Ma — 1.76. In a way similar to the case of small Ma 
and almost perpendicular propagation, the density PDF 
shown in Fig. 8, displays an excess at small density. 

At small values of Ma, the magnetic field undergoes 
mild variations. The perpendicular forcing preferentially 
excites fast waves but the slow waves that form by nonlin- 
ear interactions are the most effective at producing den- 
sity fluctuations since, the field being almost straight, only 
thermal pressure acts against compression at dominant or- 
der (see below for the effect of Alfven wave pressure). In 
contrast with the case of large Ma, a larger scatter is ob- 
served on the magnetic pressure vs. density diagram, as 
seen on Fig. 9. As seen on Fig. 10a, the correlation be- 
tween density and magnetic fleld intensity is lost and the 
density PDF is closer to a lognormal (Fig. 10b). 
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Fig. 6. Log-log scatter plot of the magnetic pressure ver- 
sus density for two simulations with forcing on the three 
velocity components and cosO = 0.7. The parameters of 
the left panel are Ms = 2.90, Ma = 0.521, P = 0.00375, 
^ = 1.20, ^ = 0.473 and a = 0.32. For the right panel 

Ms = 6.59, Ma = 2.48, = 0.141, ^ = 1.43, ^ = 9.62 
and a = 0.49. 




IE 10 0,." C E 1.0 10 12 

Fig. 7. Log-log scatter plot of the magnetic pressure ver- 
sus density for two simulations with 9 — 0. The pa- 
rameters of the left panel, where the forcing only affects 
the perpendicular velocity components, are Mg = 2.72, 

Ma = 1-76, is = 0.440, ^ = 1.04, ^ = 4.73 and 
a = 0.21. For the right panef, where the forcing is applied 
on the throe velocity components, Mg = 5.10, Ma = 2.36, 

/3 = 0.218, ^ = 1.73, ^ = 7.91 and a = 0.57. 



The behavior of the Alfven wave pressure can also be 
tested numerically. We have chosen to take (3 = (by 
turning off the thermal pressure term), in order to make 
its effect more visible. In a run with parallel propagation, 
the forcing is now applied on modes 49-51 for the per- 
pendicular velocity components and on modes 1-5 for the 
longitudinal component. Three different simulations have 
been performed, with Ma = 0.073, 0.795 and 2.09. As vis- 
ible on Figs. 11-13, while the magnetic pressure roughly 
behaves as p^/^ for Ma = 0.4, a clear tendency to ap- 
proach a polytropic behavior with 7 = 2 is observed when 
Ma increases. 



5. Conclusions 
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Fig. 8. Top: Snapshot of logj^Q p (solid) and logj^Q |6p (dot- 
ted) for the run of Fig. 7a. Bottom: Density PDF for the 
run of Fig. 7a in log-log coordinates. The dashed shows a 
lognormal fit. 

5.1. Summary and discussion 

In this paper we have investigated the dependence of the 

magnetic pressure, (and therefore, of the magnetic field 
strength) on density in fully turbulent flows, parameter- 
ized by the Alfvenic Mach number Ma (which provides 
a more direct measure of the amount of field distorsion 
than P) of the flow and the angle between the mean field 
direction and the direction of wave propagation. To do 
this, we first employed the simple- wave formalism to ob- 
tain insight by deriving relations between the fluctuations 
of the magnetic and velocity fields on one hand, and the 
density fluctuations on the other. We then presented nu- 
merical experiments confirming the expectations from the 
analysis, and extended the results. 

From the simple- wave analysis, we concluded that the 
production of density fluctuations is dominated by the 
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Fig. 9. Log-log scatter plot of the magnetic pressure ver- 
sus density for a simulation with = and forcing on 
the perpendicular velocity components. The parameters 
are = 5.34, M„ = 0.172, P = 0.00108, ^ = 1.37, 
^ = 0.168 and a = 0.56. 

slow mode at low Ma, while at large Ma both modes con- 
tribute. This result is in disagreement with the assumption 
of Lithwick & Goldreich (2001) that the fast mode can be 
neglected in accounting for the spectrum of density fluc- 
tuations because it decouples from the Alfvcn and slow 
modes due to its larger phase speed. Moreover, except for 
the case of small field fluctuations with simultaneously 
moderate density fluctuations and parallel field compo- 
nent, the magnetic pressure behaves as 

-Pmag ^ Ci - C2p (35) 

for the slow mode, and as 

Pmag ^ (36) 

for the fast mode. This different scaling of the magnetic 
pressure with density for the two modes is at the basis 
of the dccorrclation observed, as in general a given value 
of the density can be arrived at by a different kind of 
wave, and therefore have a different associated value of 
the magnetic pressure. More generally, we can say that 
the particular value of the magnetic pressure of a fluid 
parcel will not be uniquely determined by its density, but 
instead, that it will depend on the detailed history of how 
the density fluctuation was arrived at. This also implies 
that turbulent pressure cannot be simply modeled with a 
polytropic law in general, and in fact it can act as a "ran- 
dom" forcing instead of as a restoring force. In particular, 
strong density fluctuations are possible even in the pres- 
ence of a large uniform magnetic fleld, but in general there 
is no relation between Sp/p and /?. 
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Fig. 10. Top: Snapshot of log^g P (solid) and log^g |6p 
(dotted) for the run of Fig. 9. Bottom: Density PDF for 
the run of Fig. 9 in log-log coordinates. The dashed shows 
a lognormal fit. 



The angle 6 between the the mean field and the di- 
rection of wave propagation also plays an important role 
in determining the relative importance of the modes. At 
perpendicular propagation, the slow mode does not prop- 
agate, but at a,lm,ost perpendicular propagation and low 
Ma it is dominant. At intermediate angles, the distinction 
between the low- and high-Ma cases is less pronounced, 
although at large Ma a tighter correlation between density 
and magnetic field strength is observed at high densities. 
Also, in some cases (parallel propagation at moderate to 
large Ma), density peaks can even correspond to magnetic 
pressure minima. This is in fact reminiscent of pressure 
balanced structures commonly observed in the solar wind 
and possibly present in the ISM (see the simulations Mac 
Low et al 2001). 
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Fig. 11. Log- log scatter plot of the magnetic pressure ver- 
sus density for a simulation with 9 = and (3 = The 

parameters are Ma = 0.073, ^ = 3.94, ^ = 0.0564 and 
fT = 0.60. The line segment has a slope 0.5 



Fig. 13. Log-log scatter plot of the magnetic pressure ver- 
sus density for a simulation with = and /? = The 

parameters are Ma = 2.09, ^ = 3.503, ^ = 1.01 and 
(7 = 0.76. The line segment has a slope 2. 




Fig. 12. Log- log scatter plot of the magnetic pressure ver- 
sus density for a simulation with 9 = and (3 = The 

parameters are Ma = 0.795, ^ = 4.18, ^ = 0.311 and 
a = 0.63. The line segment has a slope 2. 

Our results have implications on the functional form 
of the density PDF. We found that it tends to be log- 
normal, due to the unsystematic action of the magnetic 
pressure, which allows the thermal pressure to take con- 
trol, except when the slow mode dominates density fluc- 
tuation production. In this case, the strength of the ran- 



dom forcing-like action of the magnetic pressure becomes 
density-dependent, and causes a low-density excess in the 
PDF. 

We also presented a perturbative analysis of the 
Alfven-wave pressure, which recovered all the limiting 
polytropic cases obtained by MZ95, with 7e ~ 2 at large 
Ma, 7e — 3/2 at moderate Ma, and 7e ~ 1/2 at low Ma, 
but we also mentioned a result by Passot & Gazol (in 
preparation) using this approach showing that the Alfven 
wave pressure is in general not isotropic. As pointed out 



by MZ95, the conclusion of isotropy follows from assuming 



a unique relation between the fluctuations of the velocity 
and of the magnetic field, valid for a traveling Alfven wave. 
That isotropy does not hold in general is evidence that all 
three wave modes contribute to the velocity fluctuations, 
each one with a different functional form, similarly to the 
different scalings of magnetic pressure with density for the 
slow and fast modes. 



5.2. Implications 

Our results have a number of interesting implications on 
the standard picture for the support of molecular clouds 
by magnetic fields, and on the interpretation of observa- 
tions. Observational surveys of the magnetic field strength 
in molecular clouds seem to indicate that the former is 
essentially uncorrelated from the density at low ( ^ 10^ 
cm~'^) densities, with recorrelation seemingly appearing 
at higher densities (e.g., Crutcher et al. 2002). In par- 
ticular, in molecular clouds non- detections of the Zeeman 
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effect are generally more frequent than detections, and so 
in many cases only upper limits to the field strength are 
available (e.g., Crutcher et al. 1993; Bourke et al. 2001). 
Although it has been argued that these results are consis- 
tent with the standard theory of molecular cloud support 
(Shu, Adams & Lizano 1987) through consideration of sta- 
tistical corrections to the random orientations of the field 
with respect to the line of sight (Crutcher et al. 1993; 
Crutcher et al. 2002), it is clear that they are also con- 
sistent with the lack of correlation between field strength 
and density found in the present paper and in other nu- 
merical simulations (Padoan & Nordlund 1999; Ostriker 
et al. 2001). This has been already pointed out by Padoan 
& Nordlund (1999). 

The recorrelation apparently observed at higher den- 
sities (also noticeable in the 3D simulations of Ostriker 
et al. (2001) including self-gravity) is consistent with our 
cases with intermediate values of 9 and large Ma, suggest- 
ing that the Alfvenic Mach number in molecular clouds 
is relatively large. However, it is also possible that this 
is an effect of self-gravity becoming important, which we 
have not considered here. It is also possible that even at 
high densities the magnetic field strength is highly variable 
from one core to another. More observations of the mag- 
netic field strength in high-density molecular cloud cores, 
reporting both detections and non detections, are needed. 
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